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Abstract. We determine analytically the distribution of conductances of quasi one- 
dimensional disordered electron systems, neglecting electron-electron interaction, for all strengths 
of disorder. We find that in the crossover region between the metallic and insulating regimes, 
P(g) is highly asymmetric. The average and the variance of P(g) are shown to agree with 
exact results. 
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1 Introduction 

The conductance g (in units of e 2 /h) of a mesoscopic disordered electron system is 
known to fluctuate strongly from sample to sample, or as a function of an external 
parameter such as a magnetic field or a gate voltage controlling the electron density 
[1]. In the metallic regime these fluctuations are universal and of Gaussian nature, 
i.e. the variance of g is given by a pure number independent of the specifics of the 
system, depending only on the presence (or absence) of time reversal symmetry with 
respect to orbital or spin motion (orthogonal, unitary and symplectic cases) [2]. For 
increasing disorder the fluctuations grow and are no longer universal, Gaussian and 
symmetric about the average value. When the variance becomes as large as the average 
conductance it is necessary to consider the full distribution of conductances, P(g). 
The situation is simple again in the localized regime, where P(g) is known to be a 
log-normal distribution, with variance ~< £n(l/g) > [3]. Except for numerical studies 
of finite size systems [4,5,6] little is known about the conductance distribution in the 
crossover regime. These studies suggest that the distribution is highly asymmetric 
[5], with —£nP(g) increasing like a power of g for g — > oo and like (£ng) 2 for g — > 
[6]. The shape of P{g) in the crossover regime depends on the spatial dimension, but 
appears to be compatible with one-parameter scaling, and hence universality at a true 
metal-insulator transition in d > 3. 

On the other hand, analytical results for finite systems (length L) in d = 2 + e 
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dimensions (e <C 1), where a weak disorder approximation can be applied, showed 
that the higher moments of P(g) are non-universal and diverge in the limit L — > oo 
[7]. It has been proposed, however, that these results are not incompatible with a 
universal distribution at the critical point, which was determined to be a Gaussian 
with power law tails [8]. This may seem surprising in view of the numerical results 
[4,6]. One should keep in mind, however, that in d = 2 + e dimensions the critical 
conductance at the transition is large < g > c = l/e> 1, which is deep in the metallic 
regime and hence is quite different from the critical value < g > c ^ 1 expected, e.g. in 
d = 3 dimensions. 

Here we consider the conductance distribution for the simpler case of a quasi one- 
dimensional wire of width W -C I, where I is the mean free path due to elastic 
scattering, and length L 3> L Although in this case (for orthogonal and unitary 
symmetry) all states are localized in the thermodynamic limit L — > oo, for finite 
length L the system exhibits well defined metallic and insulating regimes, and a smooth 
crossover between them. To be more precise, this is the case for a quantum wire with 
ideal leads of the same cross section, for which the perpendicular momenta at given 
energy Ep are quantized into N discrete levels, providing N channels of transport. 
The localization length £ = Nt in this case, which for N ^> 1 allows for a metallic 
regime to be realized in short wires (£ <L), whereas for long wires (L £) the system 
is of insulating character. For strictly one-dimensional weakly disordered systems the 
conductance distribution may be obtained analytically [9] , but in this case a metallic 
regime is absent. The dimensionless conductance g of a quantum wire can be expressed 
in terms of the N transmission eigenvalues Tj of the corresponding scattering problem 
as g = Yj^ =l Ti [10]. The joint probability distribution Pr({Ti}) of the Tj may be 
obtained [11] from a Fokker-Planck equation known as the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation [12], in the limit of large N. The distribution P T ({T 4 }) 
depends only on the parameter L/£. The DMPK approach has been shown to be in 
agreement with the exact formulation of the problem in terms of a supersymmetric 
nonlinear sigma model [13] . Within the latter formulation the average and the variance 
of the conductance have been calculated for all values of L/£ [14]. 

To calculate the conductance distribution P(g) from the joint distribution PT({Ti}) 
an N-fold integration is required, subject to the constraints < T, < 1 and EjTj = g, 
which has only been done in the limiting cases of L/£ <C 1 (metal) and L/£ 3> 1 (insu- 
lator). Here we describe a systematic and simple method, valid for all values of 
to obtain P(g) from Pr({Tj}) essentially analytically. We employ a generalized sad- 
dlepoint approximation, which recovers all the known results in the limiting cases, and 
provides results in the crossover regime in semiquantitative agreement with numerical 
data and with analytical results for the average and the variance of g. In particular, 
we find that P{g) for L/£ ~ 1 is given by a "one-sided" log-normal distribution for 
g < 1, with a Gaussian tail at g > 1 [15]. 

2 Generalized saddlepoint approximation 

It is useful to introduce variables Xi and Xj defined by Tj = (1 + Aj) _1 , Aj = sinh 2 Xi, 
in terms of which the conductance distribution may be represented as 
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Pig) = z\ „ 2^ 1 ^ dX ^ cx p [ - F ^ T ) 



(i) 



The "free energy" F for unitary symmetry (for orthogonal and symplectic symmetry 
the calculation is analogous) is obtained from the DMPK equation [11] 



i^ = 2£l/(A i ) + 5>(A i ,A j ) + £ T 



IT 



1,3 



+ \ 



(2) 



where u{\, Xj) is generated by the Jacobian of the integration over the transfer matrix 
elements and leads to "level repulsion". Here one may interpret V(Xi) = (£/2L)xf as 
a "one-body potential" and u(Xi,Xj) = —\{u\ +1*2), with u\(Xi,Xj) — £n | A, — Xj \ 



and 112 (Aj , Xj) = in 



as an "interaction potential" of charges at positions 



A,. In the metallic regime V(X) gives rise to a confinement of the charges in the 
regime Aj < 1 (note V(X) oc A 2 , A < 1), such that a description in terms of a charge 
density p(X) is appropriate. In the insulating regime (V(X) ~ £n 2 X, A 3> 1) the 
logarithmic repulsion between the charges dominates the potential V(X), leading to 
an exponentially large separation between the charges, of which only the one closest 
to the origin is of importance. 

To capture both aspects we keep the first eigenvalue Ai separate and represent all 
the other eigenvalues by a continuum density p(X), beginning at a lower limit A 2 > Ai. 
The free energy then takes the form 



F( P (X);X 1 ,X 2 ;t) 



= 2f™dXp(X)V tot (X) + 2V(X 1 ) 
+ rdXdX'p(X)u(X,X')p(X') 



l+Ai 



(3) 



where V tot {X) = V(X) + u(A,Ai) + 



Hoty^) — v V A 7 i" u V A i ^1) i" r+A " ^ nc integration over variables A 3 , ....Ajv m 
(1) is replaced by a functional integration D[p(X)]. The latter is done in saddlepoint 
approximation, leading to the integral equation for p(X) 



.Xr 



f 

JO 



dCMC - + «2(C + A2, C' + A 2 )]p(C' + A 2 ) = 2V tot (( + A 2 ) 



(4) 



where ( = X — A 2 has been introduced. This integral equation can be solved approx- 
imately by putting w 2 (C + A 2 , £' + A 2 ) = u 2 (C, C') + and neglecting Am in lowest 
order, which is exact in the limits (,(' ^> A 2 and (,(' <C A 2 . The leading correction 
term Au oc A 2 in the metallic regime can be treated perturbatively by replacing p in 
the integral involving Au by the saddlepoint solution for Au = 0, (the results of this 
approximation will be presented below). 

The saddlepoint density p sp (X) is found to develop negative parts for small A 2 , 
although p(X) is positive by definition. We take this as a signal that configurations 
of charges with A 2 < A c (for which p sp (X) starts to turn negative at small A) are 
unphysical and should be deleted. This is done by limiting the integration on A 2 to 
A 2 > A c + Ai. The free energy after the saddlepoint integration on p(X) is found as 



A 2 ;r) 



f 



dXV eff (X) Psp (X)+2V(X 1 ) 



l + Ai 
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where Ffg is the fluctuation part of the functional integral on p(A) and may be shown 
to depend on Ai, A2 as Fft = ^tn{\ 2 — Ai) + const. 

Since V e ff and p sp are linear functions in r, F(Ai, A 2 ;t) = F° + irF' + ^(ir) 2 F" 
is a quadratic form in r leading to a Gaussian integral over r in (1), with the result 

P(g) =y d\J d\ 2 e- s (6) 

z Jo J\l+\c 

where S = —(g — F') 2 /2F" + F°. The remaining integrals on Ai,A2 can be done 
numerically, or again in saddlepoint approximation. 



3 Results 



In the metallic regime <C 1), the relevant values of Ai and A2 are small of order 
L/t;. In the limit A x , A 2 — > we find \F"\=± and F' = £/L, whereas F° tends to a 
constant at the saddlepoint A 2 = A c + A x . Thus P(g) is given by a Gaussian centered 
at g — £/L, of variance 1/15, in agreement with known results [1]. 

We have calculated the expressions for F° and F 1 analytically up to and including 
all terms of order L/£ (the correction to F" is 0[(L/£) 2 ]). The correction to the average 
conductance in the metallic regime to this order is found as < g >= £,/L — T]L/£, with 
rj = 0.027, which compares well with the exact result [15] 77 = 1/45 ~ 0.022. There is 
no correction to the variance in order L/£, in agreement with [15]. 

In the insulating regime, L/^^> 1, the typical values of x\ and x 2 are both ^> 1. In 
fact, the requirement of positivity of the density for L/£ >, ir 2 /2 can only be satisfied 
if x 2 — > 00, independent of x\. Using the saddlepoint values of x\ = —^in^g/A) and 
F° = (f/L)jc? - x u F' = 4e- 2x i and F" = -fe" 4 * 2 -» one finds 

P(g) = ~ exp -[±(tn(g/4) + L/if] (7) 

a log-normal distribution in agreement with [3] . 

In the crossover regime on the insulating side, where £/L < 1, we make use of the 
fact that the typical values of x±,X2 are x 2 3> 1, but x\ < x 2 , otherwise arbitrary. 
We then find F° ~ (l/3)(f/L) 2 a;| - (£/L)(x 2 2 - x\) + x 2 - (l/2)£n(x 1 smh(2x 1 )): 
F' ~ cosh" 2 Xl and F" = - sinlr 2 (2x 2 )[± - ^ + sinh" 2 (2x 2 )]. 

The saddlepoint equation for X\ is given by coshx! = g^ 1 ^ 2 , which has a solution 
only for g < 1. For g > 1, instead the boundary values X\ = 0, x 2 = (2L/tt£) give the 
minimum of F. The corresponding results for P(g) are 

P(.g)-|exp- a (. 9 -l) 2 , g>\ (8) 



Z g 



arsechy/g 



1/2 



exp 



|- (arsech^/g) 



< 1 



Here a = F"(x 2 — 2L/ir£) controls the Gaussian cut-off of P(g) for g > 1. 
L/£ ^> 1 and g -1 S> 1 Eq. (9) reduces to the log-normal distribution (7). 



(9) 
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Fig. 1 Conductance distribution P(g) versus g for £/L = 0.4,0.7,1.6,2.0 (dotted, solid, 
long-dashed, short-dashed lines) 



Fig. 1 shows P{g) versus g for several values of £/L = 0.4,0.7,1.6,2.0. In Fig. 
2 the results for £/L = 0.4 and 0.7 are again shown plotted versus in(l/g) together 
with results for £/L = 0.25,0.1. In the logarithmic plot (Fig. 2) one clearly rec- 
ognizes a log- normal distribution centered at in(l/g) = L/£ — in4 and of variance 
var[£n(l/g)] = 2L/£, cut-off at ing = 0(g = 1). The abrupt qualitative change of 
the shape of the distribution P{g) in the crossover regime (L/t; <~ 1) as one goes from 
values g<ltog>lisa consequence of the small value of F". Note that | F" |<C 1 
even in the metallic regime, decreasing exponentially in the insulating regime. The 
term oc (g — F') is thus multiplied by the large number 1/ | F" |, which forces the 
saddlepoint equation F' = g, which, however, has a solution only for g < 1. This leads 
to a power law dependence of S on ing. At g > 1 the minimum of S is attained at the 
boundary of the integration regime xi — 0, where the term oc (g — F') 2 dominates, 
resulting in a Gaussian cut-off of P{g). 

This result might suggest that quite generally the statistics of the conductance in 
the crossover regime is Gaussian centered at g = 1 for g > 1 and log-normal centered 
at < ing >— L/£ for g < 1. This may be made plausible in the following way. If the 
center of the distribution is located at g ~ 1, the lowest eigenvalue Ai still must be 
dominant, meaning that A2 3> 1. The statistics of Ai is then essentially determined by 
the single particle potential V(Xi), giving rise to Gaussian statistics for Ai < 1, where 
V(Xi) oc Af and to log-normal statistics for Ai > 1, where V(Xi) oc in 2 Xi. These 
dependences carry over to the statistics of g ~ 1/(1 + Ai). 

The shape of P(g) in higher dimensions, as determined numerically [5,6] shares the 
feature of an abrupt change at g = 1 from approximately log-normal to exponential 
behavior with our results. Thus, although the DMPK approach followed here is ap- 
plicable only for quasi one-dimensional systems, the qualitative behavior found here 
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Fig. 2 Conductance distribution versus £n(l/g) for f/L 
long-dashed, short-dashed lines) 



0.7,0.4,0.25,0.1 (solid, dotted, 



may be more generally valid. 

We acknowledge useful discussions with A.D. Mirlin and D.G. Polyakov. This work has been 
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